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Abstract 

We study the anomalous dimensions and coefficient functions generated by the BFKL 
equation in 4 + 2e dimensions, by investigating both running coupling effects, and the 
inclusion of the full next-to-leading kernel. After generalising the Fourier representation of 
the solutions to this case, we analyse the /3-dependent renormalisation-group factorisation 
of anomalous dimension and coefficient contributions to the gluon density. We derive on 
this basis the normalisation factor of the Qo-scheme with respect to the MS-scheme, 
including /3-dependent corrections to it, and we outline the derivation of the full next- 
to-leading contributions. We also provide an expression for the resummed 7 qg in the 
MS-scheme which exhibits its universality and is explicit up to quadratures. 



1 Introduction 



The relationship of the BFKL equation [1] — describing small-x evolution in QCD — with the 
DGLAP equation [2] — describing Q 2 evolution — has been the subject of several investigations 
in the nineties [3-9], in the attempt of providing a unified picture of small- x physics in QCD. 
It is well known that consistency of the BFKL approach with renormalisation group (RG) 
factorisation is achieved by means of resummation formulas of large contributions at small- x to 
the quark and gluon anomalous dimensions and coefficient functions, which have been derived 
at the leading- log a; (LLx) level [3] and, in part, at the next-to-leading (NLx) one [5,6]. 

In particular, Catani and Hautmann [6] have studied the LLx BFKL equation in 4 + 2e 
dimensions with dimensional regularisation (and frozen coupling) and have derived, on this 
basis, the resummation of leading logarithms [5, 6] for the gluon coefficient function — usually 
called R(a s ) — and for the quark anomalous dimension 7 qg in the MS- and DIS-schemes. This 
calculation extends in part into the NLx level, not only because of 7 qg being next-to-leading, but 
also because of the .R-factor, which provides NLx corrections to the gluon anomalous dimension 
7 gg when running coupling effects are turned on [10]. 

It is known, on the other hand, that full NLx terms [11,12] are quite large and negative, and 
that doubly resummed approaches [13-15] are required in order to stabilise the subleading-log 
series. This raises the question of a better analysis of anomalous dimensions and coefficient 
functions at subleading level in various factorisation schemes. While the fc-factorisation schemes 
(like the Qo-scheme [16], characterised by an off-shell initial parton) have been pushed to 
the doubly- resummed level [14], the minimal subtraction one (characterised by dimensional 
regularisation, and normally used in fixed order calculations) needs yet to be extended to a full 
treatment of NLx terms. 

The purpose of the present paper is to devise a method to perform such analysis, and to 
work it out in detail in the case of the BFKL equation with running coupling, which already 
contains quite important subleading effects. The method is then extended to a full treatment 
of NLx coefficient terms, which require the e-dependence of the NLx BFKL kernel at least up 
to O (e). Unfortunately, the latter is yet to be extracted from the literature [7-9]. 

The main tool of our analysis is the generalisation to 4+2e dimensions of the 7-representation 
of the gluon density — that is, of a Fourier representation of the BFKL solution in which 7 is 
conjugate to t = log(fc 2 //i 2 ). While for e = the running-coupling equation is a differential 
equation in 7, for e 7^ it becomes a finite-difference equation which is treated in Sees. |2]and 
and in App. [X] This allows one to write the gluon density in a generalised Qo-scheme 1 as 
the product of an anomalous dimension exponential and of a fluctuation factor that we call N. 
We then show in Sec. El that the factor 9^ = R/N is due to the O (e) dependence of the LLx 
gluon anomalous dimension. This result offers an interpretation of the mismatch of R and 3Sf 
coefficients, and a hint to possible generalisations, investigated in Sees. OK 

The general case with running coupling {b > 0) is treated in Sec. E3 and is qualitatively 
similar to the frozen coupling (b = 0) case, except that the beta-function 

(3(a s ) = = ea s - ba 2 s + O (c\ 3 ) (1.1) 

has both the dimensional contribution (e-term) and the running-coupling one (6-term). We are 
thus able to explain how RG factorisation works for coefficient and anomalous dimension parts 
for b > 0, thus relating the Qo-scheme (with dimensional regularisation) and the MS-scheme in 

1 The label Qo referred originally [16] to the fact that the initial gluon, defined by fc-factorisation, was set 
off-mass-shell (fe 2 = Qq) in order to cutoff the infrared singularities. It will be shown, however, that the effective 
anomalous dimension at scale k 2 3> Q\ is independent of the cut-off procedure, whether of dimensional type or 
of off-mass-shell one. 



an unambiguous way. In particular, we confirm that the two anomalous dimensions differ by 
the quantity dlogR/dt. 

Secs.0]and|5]are devoted to the calculation of the NLx corrections to R. In Sec.|3]we treat the 
corrections due to the running coupling, which require the (known) e-dependence of the leading 
kernel eigenvalue up to O (e 2 ), and the corresponding refinement of the 7-representation and of 
the saddle-point fluctuation formalism. In Sec. 03 we calculate the remaining NLx corrections, 
due to the inclusion of the NLx kernel. Here the finite difference equation involves two steps, 
and is truncated at NLx level by an iterative procedure. The final result involves the O (e) 
corrections to the NLx kernel eigenvalue, which are not yet explicitly known. 

In Sec. E]we reconsider the calculation of 7 qg in the MS-scheme. Contrary to the case of the 
DIS-scheme, in which the coefficient function C qg is set to zero, no closed form resummation 

formula is yet available for 7 q ^ IS ' ) . Catani and Hautmann are able to provide a recursive method 
in order to disentangle C qg from 7 qg , so as to calculate a number of terms of their expansion 
in a s /uj, oj being the Mellin moment conjugate to x, a result further improved in the later 
literature [17]. Here we use the 7-representation of the gluon density in order to do the same 
calculation: this allows us to provide an explicit resummation formula which exhibits the 
universality of 7 qg and involves only quadratures of functions which are known in principle. 
However, the latter calculation requires an all-order e-expansion of the coefficient function of 
the gluon. An interesting byproduct of the universality analysis is the proof (App. [Bj) that the 
off-shell P qg splitting function introduced in [6] is indeed process-independent. 

In the final Sec. [7] we summarise our results and we discuss their relevance for the improved 
resummations of splitting functions [14,15]. 

2 7-representation and _R-factor in the LLcc case 

Investigating the BFKL equation in 4 + 2e dimensions is of importance on its own, because - 
in addition to the ultraviolet (UV) role of dimensional regular isat ion — a positive e parameter 
allows to regularise the mass-shell singularities and, in some regime (see Sec. EJ) it avoids the 
Landau pole. For this reason we shall consider e > as an infrared cutoff, in alternative to 
setting off-mass-shell the initial condition for partonic evolution, as done in the Qo-scheme [16]. 
Our purpose is, at large, to understand the relationship of such two kinds of schemes, whether 
or not a minimal subtraction is assumed in the partonic densities. 

Since the gauge coupling g is dimensionful in 4 + 2e dimensions, we shall introduce, as usual, 
the renormalisation scale fi, the dimensionless coupling in the MS-scheme 



and the parameter t = log(fc 2 //i 2 ) in terms of which the running coupling is a s (t) = a s e £t . 
Therefore, the LLx equation shows a running coupling with infrared free evolution correspond- 
ing to the dimensional contribution of the beta function (jl.lj) in the b — > limit. Since a s (t) — > 
for t — > —00, we can set on-shell (feo = 0) the initial condition for the gluon Green's function, 
so that e acquires the role of infrared cutoff. The ensuing solution of the LLx BFKL equation 
is naturally expressed as a power series in d s (t), which was extensively studied in [6]. 

Our purpose here is to recast the solution for the gluon density mentioned above in the 
7-representation form — 7 being conjugate to t — so as to be able to describe in a simpler way 
its anomalous dimension behaviour. The BFKL equation for the unintegrated gluon density JF 




(2.1) 



in 4 + 2e dimensions reads (the cu-dependence of T is understood) 
F E {k) = 5^\k) + i J ^JL X £ (k, k')F E (k>) 
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X £ (k,k')=g 2 N i 
and its power series solution is 
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is the LLx "characteristic function" of % £ in 4 + 2e dimensions, defined by 
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The series (|2.4|) is well behaved in the infrared (t — > — oo) apart from the \je poles, but is 
not really suitable in the ultraviolet (t — > +oo), where the variable a s (t)/u grows, and at most 
a finite radius of convergence is expected. We thus look for a solution in the 7-representation 
form 



Fe(k) = I e*/.(7) 



»7=c 



2vri 



(2.7) 



where the Fourier variable 7 has the interpretation of a continuum in which 7 = ne is the 
lattice counterpart. We expect Eq. (|2.7|) to be best suited for anomalous dimension properties 
in which n — > 00 and £ — > with 7 = ne kept fixed. 

Let us start noticing that the "ansatz" ([2.7)1 . when replaced in the BFKL equation with a 
general initial condition f E (j), leads to a finite difference equation of the form 



/ £ (7 + e) = /i 0) (7 + ^) + -X £ (7)/ £ (7)- 



In the following we shall often consider the case 

/f ( 7 + e ) = /(o)( 7 )^ 
which corresponds to the initial condition 

5f\k) 
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i.e. to the one expected from Eq. ()2.4j) with a cutoff at k 2 = Ql = u 2 e~ T . In this way we can 
study the role of the cutoff, by solving (|2.8J) with T fixed, and performing the limits T — ► +00, 
£ —>■ either in this order, or in reverse order. 

We shall solve Eq. ()2.8j) in two steps, starting from the homogeneous equation 



h e (j + e) 
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Xe (7)/ ie (7)=e^«/ i£ ( 7 ) 
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where we take the ansatz 



h^) = exp 1 1 7 £ £ ( 7 ')^ j = exp (S £ ( 7 )} (2.12) 
so that Eq. (j2.11|) implies 

5 e ( 7 + e) - 5 £ ( 7 ) = f +£ £ e ( 7 0— = ^(7) • (2.13) 

We show in App. IA. II that the solution for the "Lagrangian" C £ (j) is expressed, under appro- 
priate conditions, in terms of the Bernoulli numbers B n as follows 

00 „ 



n=0 



where L^ n > is the n-th derivative of L with respect to 7, and the generating function is 

00 „ 

y^ z » = _f_ (2.15) 
^ n! e 2 - 1 v 7 

n=0 

so that 

B = l, ^ = -1 52 = 1, ,63 = 0, B 4 =~---. (2.16) 
2 o d(J 

Correspondingly the "action" takes the form 

5.(7) = J / 7 - ^(7) + ^4(7) + O {e 2 ) , (2.17) 

where L £ can be further expanded in e as follows 

£,<?) . log (*x.( 7 )) = log I4M) + ( £ 3) (2 . 18) 

V ^ 7 V ^ ' Xo(t) VXo( 7 ) 2x0(7)/ 

and the BFKL eigenvalue function % £ = Xo+ £ Xi+ £2 X2+C (^ 3 ) is given, according to Eq. (|2.5j) . 
by 

Xo(7) = 2V(l)-^(7)-^(l-7) (2-19) 
Xi(7) = 5 [%o(7) + 2^(1) - ^(7) - il/ (1 ~ 7)] (2.20) 
X 2 (7) = Xo(7) [Xi(7) - ^'(1)] - 1^(7) + ~ [8^(1) - V"(7) - <(1 - 7)] ■ (2.21) 

Therefore, the solution of the homogeneous equation in t-space (Eq. (12. 2j) without delta term) 
takes the form 

where we have truncated the e-expansion of the exponent up to the finite terms. The choice of 
the lower bound of the 7' integrals is of no concern, since the normalisation of the solution of 
the homogeneous equation is arbitrary. 



The inhomogeneous equation (|2.8|) has, on the other hand, the iterative solution 

fed) = / (0) (7 -e) + / (0) (7 - 2e)-x £ (7 - e) + / (0) (7 - 3e)^x £ ( 7 ~ 2e)-X e ( 7 -£) + ••• 

o> a; a; 

00 - n— 1 n— 1 00 h ( \ 

=E/ ro (7-«)(^) n^-")'i^-») u 7 rL) ' (2 - 23) 

n=l m=l n=l eV ' y 

which is closely related to the power series solution (|2.4J) . as analysed in more detail in App. IA.2l 
Here we just note the e = limit of Eq. ()2.23|) at fixed T. Due to the fact that ^(7) satisfies 
Eq. (|2.11|) . the ratios of h's has a non trivial e — > limit, and we obtain 

/o(7) = / (0) (7)Ee nLoW = 1 /( ! ( v 7) r , , (2.24) 

as expected for the solution of the BFKL equation with frozen coupling and with cutoff. We 
get in fact, up to higher twist corrections, 

MQ= n , (7-2) 0(^ 2 -A- T ), (2.25) 
-7oXo(7o) \fi v 

where 7o(a s /c<j) is the LLx anomalous dimension, defined by 



-Xo(7o), 7o - =- + - • (2.26) 

UJ \UJ / UJ \ UJ J 

The result (I2.25|) shows the customary infrared singular T-dependence of the Q - sc heme, de- 
termined by 7 . 

Our emphasis in this paper is, however, in the limit of (|2.23|) for e <C 1 and T — > +00, where 
we obtain its continuum counterpart which roughly replaces the sum in the r.h.s. of (|2.23|) with 
an integral, apart from some e-dependent corrections. In order to specify such corrections, we 
rewrite Eq. ()2.8|) in terms of p £ {l) = fe(j)/h E (j) in the simpler form 

^ +E) -^ ) = ^^ = w e ~ i " (7, • (2 ' 27) 

which is of the type (j2.13j) and is solved in terms of Bernoulli numbers as in (|2.14)1 . The r.h.s. 
of ()2.27|) has, however, 1/e singularities in the exponent (cf. Eqs. (j2.12|) and (j2.17|0 . and its 
derivatives generate eventually a non trivial correction factor (see App. IA.3|) having a finite 
e — > limit, as follows: 

nM /(0)(V) L ^- ST x n I AMI f9 98^ 

where we have kept terms ~ eT, that we consider a number of order unity. 
Finally, by replacing Eq. ()2.28|) into (|2.7|) we get the solution 

~ r d 7 r' , / exp l 7t + « # ^(^)d, + iT - 1 Lo(*)d* + /; ggd,} 
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where we have truncated the e-expansion to the finite terms. This solution has the form of a 
double 7-representation, similarly to the customary case b 7^ 0, e = 0, in which the e-parameter 
replaces b by providing an infrared free running coupling. 



Our goal is now to investigate the T — > +00 limit of Eq. ()2.29|) by removing the artificial 
infrared cutoff Qq = fi 2 e~ T . In this limit the leading 7'-dependent phase is generally large and 
given by 

E{i) = 7T - - T L (z)dz , (2.30) 
£ Jo 

so that the 7'-integral is dominated by the saddle point 7': 



which implies 



E\i') = T — ~L (j') = , (2.31) 

1 = — -Xo 70 • 2.32 

OJ 

Therefore, 7' ~ a s (— T)/u; — > for T — ► +00 2 and, by taking into account also the 7'- 
fluctuations 

1 (233) 



7 " V^W) V -^o(y) ' 

all T-dependent factors in ()2.29|) cancel out, yielding 

Mk) T ^° I d J_ J__ e W \ 1 t+ 1 - / 7 L (y)dy+ r^) d y + o( £ )| , (2.34) 

J V2^i Vxo(7) I £ Jo Jo Xo(y) J 

which is just a solution (|2.22j) of the homogeneous equation with an appropriate normalisation. 

Eq. (|2.34j) is the main result of this section and shows the mechanism by which the solution 
becomes independent of the details of the initial condition in the T — > +00 (Qq — > 0) limit. In 
fact, due to the infinite evolution path from —00 to t, the shape of the solution reduces to the 
one of the homogeneous equation, and the initial condition determines only the normalisation. 
This already suggests the factorisation of the 1/e singularities which have replaced the In- 
dependence in Eq. (I2.34|) . In order to prove such factorisation we resort again to the saddle 
point method in order to evaluate the e —>■ behaviour of ()2.34|) . The stationarity condition is 
now 

et + L (7) = i = ^l Xo (7), (2.35) 

with a stable fluctuation along the real axis for x'o(7) < 0; thus yielding the result (7 = 7^ = 

7o(a 8 (t)/w)) 

= i=exp( 7 ^+i ^ L,{iW+ r^iMdylx [l + 0(e)]. (2.36) 

V-x'o(it) I e Jo Jo xo(y) j - 

Finally, we calculate the integrated gluon density in the form 
g £ (t) = J d 2+2e fc' F £ {k') 0(fc 2 - k' 2 ) 

= J-— e X pUt + - L (y)dy+ r^dyjx [l + 0{e)] } (2.37) 
W-Xo(7t) I £ A) Jo Xo(7') J 

where we have performed an integration by parts of in (|2.36|) . 

2 The solution with 7' ~ 1 is relevant in the infrared region k 2 <C Qq (t <C — T), where it behaves as 



Some remarks are in order. Firstly, the saddle point condition ()2.35|) provides a real positive 
anomalous dimension only in the perturbative regime in which a s (t) < uj/xq{^)i value at which 
the LLx anomalous dimension has the well known singularity at 7 = ~. Since, however, 
a s (t) increases with t for e > 0, for large t the solutions to (|2.35j) will become necessary 
complex conjugate with 9^7 = |, and g e (t) will oscillate asymptotically This shows how the 
perturbative series (|2.4jl behaves beyond its convergence radius, so that T{k) stays marginally 
square-integrable, 3 as required by the 7-representation. 

Secondly, Eq. (J2.37|) shows, in the perturbative regime, a fluctuation factor N and an anoma- 
lous dimension exponential, in the form 
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is the O (e) correction to the BFKL anomalous dimension. This interpretation follows simply 
by identifying the exponents in (J2.38|) and (J2.37j) . which have the same t-derivative (by the 
saddle-point condition (|2.35|1 ) and the same value at t — > —00 (jt = 0). Note that the 1/e 
singularity of the exponent is a consequence of the boundary at t — > — 00 in ()2.38|) . because of 
the Jacobians 

d« s Xo(7t) j- 
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which relate the various forms of the 1/e exponent: 
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We are now in a position to factorise the t-dependence from the 1/e singularities in Eq. (|2.38|) . 
By simply subdividing the r integration interval into ] — 00, 0] U [0, t], we obtain 4 
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where we have defined the reduced coefficient factor 
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which arises because of the anomalous dimension correction e^i cancelling the 1/e singularity 
of the Jacobian (j2.41|) according to the identity 
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3 That is, .F(fc) ~ (A; 2 ) 1 / 2 x oscillating function. 

4 Splitting the integration interval at r = corresponds to choose the factorisation scale [if = fj,, i.e., = 0. 



The expression flCTl) agrees, by Eq. with Eqs. (3.17,B.18) of Ref . [6]. Finally, the finite 

t-evolution in Eq. (|2.43jl becomes simply 




(2.46) 



in the e — > limit, thus recovering the result [6] 

This derivation shows the mechanism by which the R coefficient factor is obtained as the 
product of N (arising from the 7-fluctuations) and 9^ (arising from the ^-dependence of the 
BFKL anomalous dimension) . In the following we shall generalise this mechanism to b > and 
to further terms in the e-expansion. 



3 The running coupling equation and its factorisation 
properties 

Our first purpose is to generalise to 4+2e dimensions the BFKL equation with running coupling. 
Because of the e-dependence of the /3-function in Eq. 1)1. ip . the running coupling a s (t) has the 
form 

1 b _ £t [l V 
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where b = 11/12 in the Nf = limit. Due to the UV fixed point of Eq. (jl.ip at a s = e/b, 
Eq. (J3.1j) shows two distinct regimes, according to whether (i) a s < e/b or (ii) a 8 > e/b. In 
the regime (i) ct s (t) runs monotonically from a s = to a s = e/b for —00 < t < +00, while 
in the regime (ii) a s (t), starting from e/b in the UV limit, goes through the Landau pole at 
£a = log(l — e/ba s )/e < 0, and reaches a s = e/b from below in the IR limit. 

Since the LLx kernel % £ (t,t')dt' scales like e £t , an equation which realises such a s (t)- 
evolution (and regimes) is obtained by setting 
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It is soon realised that JF e 6 (fe), given by (cf. Eq. ()A.5 
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has a well-defined iterative solution in the regime (i) in which < a s (t) < e/b. In fact, the 
fc-integrations are convergent in the IR because of a s (t) ~ e et for t — > —00, and everywhere else 
because a s (t) is bounded. In this sense, b, e > in the regime (i) act as regulators of both the 
IR and UV regions, and of the Landau pole, so that Eq. ()3.2j) is meaningful without any cutoffs 
- a somewhat unique case in BFKL theory. In the b — > limit Eq. (J3.2j) reduces to Eq. ()2.2p . 
but become less convergent in the UV region: as noticed in App. IA.21 only a finite number of 
iterations n < 1/e is actually possible for b 7^ 0, for the UV integrations to be convergent. 



Therefore we shall solve Eq. (|3.2|) in the regime (i), where it is nicely convergent, and we 
shall study the factorisation of 1/e singularities by letting b, e — > with e/b > a s kept fixed. 
Eventually, we are instead interested in the physical limit of e — > with a s and b kept fixed, 
but we shall perform it only at the end, after factorisation of the 1/e poles. 

Let us introduce the 7-representation ()2.7|) into Eq. ()3.3|) after multiplying it by the in- 
dependent denominator. After simple algebra we get the equation 
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where we have again introduced the cutoff t > —T in the initial condition (cf. Eq. (|2.10J) ). in 
order to better control the e — > limit. Eq. (|3.4jl has the form of Eq. (J2.8)) and admits a similar 
iterative solution, which is obtained from Eq. ()2.23j) by the replacements 
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(3.5b) 
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In the limit of b, e — > with b/e and the cutoff T kept fixed we get identically from Eqs. ()2.23j) 
and ()3.5|) the frozen coupling solution ()2.24|) and ()2.25|) . as expected. 

We are more interested, however, in the continuum limit of the solution of Eq. ()3.4j) in the 
regime (i) and for T — ► +00. By formal manipulations similar to the 6 = case we obtain 
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which is again a solution of the homogeneous equation, where 
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and the exponent has been expanded in e (at ba s /e fixed) up to the finite terms. Once again, the 
factorisation of the 1/e poles in ()3.6|) is investigated, for e — > 0, by the saddle point condition 
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where, due to the 6-dependence of Eq. (|3.7j) . a 8 (t) has the expected form (|3.1j) . It follows that 
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(3.10) 



The final expression in Eq. ()3.10|) is formally identical to Eq. (|2.38|) and is expected to have 
analogous factorisation properties. However, due to the different form of a s {t), the Jacobians 
induced by Eq. (|3.8|) are different: 



dt 



da. 



Xo(7t) 



-dj, 



t , 



ba s ) e [ X o(7t) " t] 
and this explains the different e-dependence of the exponents at fixed value of a s (t): 



dr 7q 
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(3.11) 
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(3.13) 



/ 1) ill. — ba ) 

The factorisation of the t-dependence in Eq. (I3.10|) from the 1/e singularities is now per- 
formed as in Sec. 121 by subdividing the r-integration into the ] — oo, tf] and [tf,t] intervals, 
where tf = log(fij/fi 2 ) defines the factorisation scale. The two intervals are treated differently: 
in the finite UV interval [tf,t] we can freely go to the e = limit at fixed value of a s (tf), thus 
recovering the a s (t) —> a B (tf)/[l + ba B (tf)(t — tf)] limit and the normal t-dependence of UV free 
QCD in 4 dimensions. In the remaining infinite IR interval we have to factorise the 1/e poles 
by expanding in the ba s /e parameter, as normally done in fixed order perturbation theory. 

By looking at Eqs. (|3.12|) and (|3.13j) (with t replaced by tf) we realise the following: 
Eq. ()3.12|) exponentiates 1/e poles and higher order ones coming from the 6a s /£-expansion, 
and has to be factorised in full in the MS-scheme; it evolves in tf according to the LLx anoma- 
lous dimension. On the other hand, the term (|3.13j) - - due to the O (e) correction to the 
BFKL anomalous dimension — reduces, in the b — > limit, to the ^Pi{a. s {t f) / ui) factor found in 
Sec. 121 which yields RfN, i.e., the normalisation mismatch of the frozen coupling evolution with 
respect to the MS density. Furthermore, in the fcs/e-expansion, Eq. ()3.13j) exponentiates 1/e 
and higher order poles which should be factored out in the MS density, and therefore contribute 
to the MS anomalous dimension. 

In order to perform collectively such separation, we rewrite the r.h.s. of Eq. ()3.13|) (with 
t = tf) in the form 



^ {t f ] da (a 

— 7i-l 

a \uj 



a{e 



da 



— ba 7i ( - 
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(3.14) 



where the first term provides logD^as^/)/^), and the second one is the anomalous dimen- 
sion contribution to be factored out in the MS density. Finally, by replacing Eq. ()3.14|) into 
Eq. ijTTTty . we find 
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where 



^ MS )(t / ) = exp 




7o( — J + 6o7i( — 



(3.16) 



incorporates all the 1/e poles, in a formal 6ft s /e-expansion. 

Eqs. f!3.15|) and ()3.16|) are the main results of this section. They confirm the relation between 
our generalised Qo- scneme (with dimensional regularisation) and the MS-scheme, namely 



g( Qo \t) = g s (t) — R 



a s (t) 



UJ 



g {MS) (t) , 



(3.17) 



and they prove the following resummation formulas 



d_ n (MS) 



y» Wfliw) .^-^) + * w g5 (3 , 9 ) 

=70( M) )+6a , (t)71( ^) 

where N, *H, 7 and 71 are defined in Eqs. lj£33)l . (|2~H|l . (jQgj) and (l2~471|) respectively. 

Let us remark that the results (|3.18H3.iTH) are not surprising, in view of the well studied 
relationship of the N and R factors already performed in the literature in a somewhat different 
context [10]. However they offer some important insights into the explicit RG factorisation 
obtained here: Eq. (j3.18|) shows that the evolution of the gluon density in a fc-factorisation 
scheme is independent of whether the BFKL equation is regularised by a cutoff (Qo-scheme 
in strict sense) or by a positive e (as done at present). This is a consequence of the RG 
factorisation just proved. Eq. (j3.19|) shows that the MS evolution is the one expected because 
of the R factor in Eq. (J3.17|) . However, this result is obtained by the explicit, 6-dependent 
factorisation of the 1/e poles connected with 7! in Eq. (J3.16j) . Furthermore the t ^-dependence 
in Eq. ()3.15|) cancels out because of the 71 evolution in g( MS ) cancelling the ^-dependence of 
9^. All together, such terms build up the O (e) correction to the BFKL anomalous dimension, 
as done in Eq. (|3.14j) . which precisely vanishes in the e — > limit. 

Our next goal is to generalise the above factorisation procedure to subleading terms, in 
particular to the coefficient terms at NLx level and to the corresponding NNLx anomalous 
dimension terms. 



4 Factorisation at subleading level: 6-dependent correc- 
tions 

We have just realised that the normalisation factor arises from the O (e) correction to the 
BFKL anomalous dimension (|3.13|) after the (minimal) subtraction of a 6-dependent contribu- 
tion to the MS anomalous dimension (second term of (j3.14j) ). This suggests that, in order to 
compute corrections of relative order ba s to and the corresponding ones to 7^ MS ) we have 
to calculate O (e 2 ) corrections to the BFKL anomalous dimension and, more generally, to the 
exponent of the solution ()3.6|) of the homogeneous equation. We thus restart our analysis from 
the "action" (j2.17J) and its 6/e-dependent counterpart, and we rewrite the solution ()3.6|) of the 
homogeneous equation in the form 

*U*)= / %= 1 1 =exp( 7 t+- f 1 Ls,HW ' + ^L'eM + '(e 2 )| > 

J V27r£ \JXe(l)- b -f 1 £J ° 12 J 

where we have kept the ^-dependence of L and x> an d the O (e) correction to the "action" 
(which is of relative order O {e 2 )). 

We then look at the factorisation properties of Eq. (j4.1J) by expanding around the e- 
dependent saddle point 7 = 7 e : 

et + L £ , b {%) =0 ^ i = ^l Xe (7 £ ), (4.2) 



which, therefore, defines the e-dependent BFKL anomalous dimension 
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where 71 was given before in Eq. (|2.40j) . and we have 



72 — 
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(4.4) 



7=70 (ff) 



Furthermore, we have to compute the 7-integral by expanding the action and the factor 

[Xe(t) — ~rl around the saddle point. Since the 7-fluctuations are governed by the width 

a 7 = [— e/L' Eb (%)] which is of O (y/e) while the "action" is 0(l/e), we need to expand, 
it turns out, up to 6-th order in 7 — 7 in order to reach all O (e) terms. This calculation is 
performed in App. [Oand, when replaced into Eq. (j4.1|) . provides the following result for the 
gluon density: 
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Here we have used the identity 
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the Jacobians 
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and the O (e) correction to the action 
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(4.8) 



where the terms in square brackets are precisely the fluctuations calculated in App. IC.ll 

Let us now look at the factorisation properties of our result in Eq. (J4.5)) in the limit of 
e, ba s /e — > 0. The anomalous dimension exponential is factorised as usual, and its infrared part 
at the factorisation scale tf is given by 
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We have now the additional O (e 2 ) term 72, which however still builds 1/e singularities because 
of the 6a s /e-expansion of the /3-function in the denominator. The decomposition into coefficient 
and MS-anomalous dimension parts is simply done by writing e 2 = (ba s ) 2 + (e — ba s )(ba s + e), 
so that 
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(4.10) 

where the first term in the r.h.s. is the minimal subtraction part of the 72 anomalous dimension, 
the second one is the O (ba s ) correction to log 91 



We have also to look for 1/e singularities possibly arising because of the 6a s /e-expansion in 
the contribution (|4.8|) . Here the situation is not so clear a priori, given the fact that 
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(4.11) 



has a non-trivial ba s /e expansion. However, after some algebra, we find a remarkable cancella- 
tion leading to the result (see App. IC.2J) 
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which is only linear in b/e. It follows that S\ has no minimal subtraction terms, and only 
contributes to the renormalisation of the N-factor, as follows: 
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(4.13) 



We can thus summarise our results on the NLx corrections to the coefficient factors of 
Eq. (j3.15J) (and on the NNLx ones to the MS anomalous dimension). The 91 factor takes the 
form 
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and the MS anomalous dimension becomes 
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d log a s 



Due to the form of (|4.15|) and (|4.14|) . the t/-dependence of the factorisation formula (|3.15|) 
cancels out as it should, and as is expected on the basis of their common origin, Eq. ()4.9|) . 
We also expect this mechanism to hold true to all orders in ba s , so that the MS anomalous 
dimension and coefficient can be inferred from the corresponding e-expansion of the BFKL 
anomalous dimension (J4.3|) . 

On the other hand, the fluctuation factor 74 takes finite NLx corrections from Eq. (J4.13|) . 
and becomes 
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(4.16) 



Such corrections are identical to those found from the normal 7-representation in 4 dimensions. 
In this case also we expect that higher orders in the ^-expansion ()2.17|) and the corresponding 
fluctuations will combine so as to provide further finite subleading corrections to Eq. ()4.16|) . 
We have no formal proof of this expectation, but we notice that higher order fluctuations still 
involve the scale t and the coupling a s (t) and are thus independent of the IR behaviour at 
t — > — 00. It is then natural to believe that they can be computed directly in the e = limit 
from the normal 7-representation in 4 dimensions, as explicitly proved above for the next-to- 
leading terms. 



5 On the treatment of the full NLcc corrections to the 
£K factor 



Having calculated the running coupling corrections to fR, the problem arises of including all 
NLx contributions, as embodied in the next-to-leading BFKL kernel. The structure and the 
explicit expression of such a kernel with dimensional regular isat ion were given in several papers 
for the gluon [7] and the quark [8,9] parts, and summarised in Refs. [10-12]. It has the form 



- (1 - e rf ) % £ 



e i 



(5.1) 



where we have singled out the running coupling part (proportional to the leading kernel % e ) 
which arises from the O (a 2 ) expansion of the running coupling equation (|3.2|) . The remaining 
kernel which scales as e 2ei , is the properly called NLx kernel, whose eigenvalues have 

been worked out in the literature in the e — > limit, and are here required up to O (e) for the 
complete calculation of NLx corrections to 9t. 

We shall thus generalise Eq. (J3.2j) to include the NLx kernel in the form 



uo 1 + ba s ^ 



{ 2+2e k > 



(27T) 



2+2e 



5 (2+2e \k) + 



^(i) 



Fs(k) . 



(5.2) 



where NNLx terms and further subleading ones have been freely added so as to reproduce the 
resummed a s (t) evolution. 5 It is now straightforward to go to 7-space and to obtain a modified 
form of Eq. (|3.4j) . By introducing the NLx "characteristic function" Xe °f in 4 + 2e: 
dimensions 

L ^(fc^o^r 1 - 2 ^^^^^-, (5.3) 



(27T) 2+2£ £ 

the corresponding homogeneous equation reads 



/ £ ( 7 + e) - ^ {W + s)- / £ ( 7 )] = a s ^M/ £ ( 7 ) + ^xP(7)/,(7 - e) 

£ UJ UJ 



(5.4) 



and contains, therefore, two finite difference steps, due to the different scaling properties of 
the leading vs. next-to-leading kernels. However, to NLx accuracy, we can replace the leading 
order equation in the last term, and we obtain 



Xeil) _ b , X^ 1) (7) 
w £ X«(7~e) 
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/ e ( 7 ) = exp[Lf 6 ( 7 )]/ £ ( 7 ) 



(5.5) 



where the next-to-leading term is now suppressed by a factor of uj with respect to the leading 
one, much in the same spirit as the u;-expansion [18]. 

We can thus solve Eq. ()5.5j) by the same method used before to get Eq. (|4.1|L Correspond- 
ingly, the saddle point at 7 = %(a s (t);uj) such that 







a s (t) 



X { e\le) 
Xe(le-e) 



(5.6) 



5 Here we keep terms which are of relative order a s e n (n > 0) with respect to the leading ones, and we drop 
terms of relative order ba^ or higher, at fixed values of e. 



admits the solution 

, E = 7 (o,(M)) +a , (t)7 a,(^)) +0(a;) , (5 . 7) 

where defines the LLx BFKL anomalous dimension of Eq. ()4.3|) and 7 e is its NLx cor- 
rection, obtained by truncating the expansion of the saddle point position to relative order 
O (<*.(«)): 

7 a, (M>) = _ x?'(ff ) x.^) = M + £7 a, ( gg) ) + H . ( , 8) 

This expression reduces to the customary one [11,12] in the e = limit, but has O (e) corrections 
coming from the corresponding ones of XeP — J e ^ to De extracted from the various papers in 
the literature [7-9]. 

Finally, we expand the e-dependence of the BFKL anomalous dimension, including NLx 
terms, as follows: 

«> - & (!) + ^ (!) + < (I) + * N" (!) - ^ (!)] («) 

and we replace it into the analogue of Eq. (|4.5jl . The corresponding anomalous dimension 
exponential factorises in the form 

exp I ^ ^ 7e("^) |> =^(«s(t);^)^ (MS) W , (5.10) 

where 

« (a> ( t ); U )=«pjy;^ [ 7 r(ewew 3] [ ( b.id 

and 

jlMS,(t) = exp {i w^)^'^ ■ (5 ' 12) 

It follows that the MS gluon anomalous dimension is 

= 7o + «s7o + & «s^j — 5.13 

a log a s 

and has therefore sub leading terms determined by the e-dependence of 7 £ at e = bd s , which are 
related in the expected way to the coefficient 9^. 

On the other hand, the gluon density in the Qo-scheme contains an additional fluctuation 
factor N(a s ; uj) whose calculation at full NLx level proceeds along the lines of Sec. Eland is not 
explicitly done here. The corresponding anomalous dimension is 

7 wo; a s ; u) = % J — J + a s % — I - ba s — — 5.14 

\u J \uj J a log a s 

and can be calculated directly at e — 0. 

It should be remarked that the NNLx contributions to the anomalous dimension in Eqs. ()5.13|) 
and (|5.14jl coming from the NLx corrections to D\ and 3\f are of course mixed with dynamical 
contributions whose calculation has not yet been attempted in the literature. Nevertheless, we 
do compute here — once 7 £ is known — the full NLx contributions to fH, which therefore 
explain the difference between MS- and Qo-scheme anomalous dimensions at NNLx accuracy. 



6 Universality of 7q^ S ^ and its resummation formula 

So far, we have considered the gluon channel only, and discussed NLx corrections to the normal- 
isation change of the gluon density from the MS-scheme to the Qo-scheme. However physical 
probes are coupled to quarks, which enter the BFKL framework through a fc-factorisation ker- 
nel He P \ defining the measuring process at hand (labelled by the superscript p), as shown in 
Fig.d 



k 










k' 







F 



Figure 1: The leading small- x contribution to the quark density. 

In a physical scheme, like the Qo-scheme, we can just define the quark density by the 
action of 

= <t) J d^k' ffW (^)^(fe') , (6.1) 

(a s (t) = a s e et , see Eq. (12.1(1 ) and it is then pretty easy to find a 4-dimensional NLx resummation 
formula for 7^ , as first shown by Catani and Hautmann [6] in the DIS-scheme (g( DIS ) = F2). 

On the other hand, in the MS-scheme one has to disentangle the coefficient part of Eq. (j6.1() 
from the anomalous dimension part. The latter contains, by definition, the minimal 1/e singu- 
larities which are not exponentiated in the gluon density in the following collinear factorisation 
formula 

4 P) (t) - Cg (a s (t)) g r\t) = gp)(t) = f dr 7 £=> (a s (r))^(r) , (6.2) 

J —00 

where we have set Cqq = 1 and 7 qq ' IS ' ) = at NLx anomalous dimension accuracy and we have 
omitted, for notational simplicity, the 10- variable. Note that C$ is process- and e-dependent, 

while 7 qg is universal and e-independent. 
q (p) i s 

obtained in the 7-representation framework by replacing in Eq. (jfi.lj) the characteristic 
function of the kernel 



k' 2 \k 2 J £ \k' 2 J 7(7 + e) 

where the two poles are expected because of the intermediate /-integration in (k > I > k' 
roughly) and the ^-dimension of the kernel integration providing also a running a s (t) = a s ■ 

c\k' 2 fk' 2 \ 1 f r\ 2+2e l (k 2 Y 



k' 2 \k 2 J J I 2 v ' v ' 7(7 + e) 



The pole I/7 in ()6.3|) is incorporated into the integrated gluon density g e (l) = fe(l)/l to yield 



qj?>(t) = a B (t) / d 7 e^ £ ( 7 ) 
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One could try a saddle point evaluation of q around 7o(o; s (t)/a;)), but it is difficult to get 
the right accuracy, and the fluctuation algorithm is cumbersome. We prefer to extract 7 qg 
directly from the expression 



#°(*) ee ±qM(t) = a s (t) / d 7 e^5£?)( 7 )ft(7) 
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(6.7a) 
(6.7b) 



having used in the first line the equality 

d 



dt 



[a s (t)e*] = a s (t)e*( 7 + e) 



(6.7c) 



(6.8) 



It is not obvious how to extract from Eq. ()6.7|) resummation formulas for both Cqg and 7 qg 
in terms of the fc-factorisation integral of qf \ For instance, at e = we get the well known [6] 
fc-factorisation result 
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which determines 7 qg in the p-scheme, but not Cqg and 7 qg separately. In Ref . [6] , a double 
expansion of (|6.7bj) in both a s (t) and e is used to get half a dozen terms in the a s (t)-expansion 
of 7 qg , a result further improved in [17]. 

Here we want to derive 7 qg 1S (a s (i)) to all orders in a s (t)/u>, by exploiting the property of 
the coefficient part in (|6.7b|) of being a total t-derivative of the product of g £ and a function 
which is perturbative in both a s (t) and e. We start noticing that "K^ can be expanded around 
7 = — e in the form 



ftW( 7 ) =ftW(-e)+0(7 + £) 
= (e) + O (7 + e) , 



(6.10) 



where 0K(e) is now process-independent, being the residue of the function (J6.1|) at the collinear 
pole 7 = — e. In App.|B]we evaluate fH^ > ^( 7 ) in various cases and we show that in all cases 0H,(e) 
is process-independent and can be written as the product of a rational and a transcendental 
part 
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This form basically follows from the off-shell generalisation of the g — > qq DGLAP splitting 
function introduced in [6], and proved in App. IB. 31 to have a universal off-shell dependence 
induced by fc-factorisation. 

On the other hand, by Eq. (|6.8|) terms of order (7 + e) n (n > 1) are easily seen to be total 
t-derivatives, so that we can rewrite Eq. ()6.7|) in the form 

M&a^R^^ie^g^it) = 7q f } (a s (t)) gf® (t) + total t-derivative . (6.12) 

This result shows that 7 qg^ S ' 1 is universal, being dependent on the universal functions ^K(e) and 
R. Furthermore, it suggests how to extract the anomalous dimension contributions from the 
power series in e in the l.h.s., i.e., by subtracting a total t-derivative or, in other words, by 
doing an "integration by parts" . By writing the expansions 

^ ran(£ ^M> )£ ) = J^i-e^f^) (6.13) 

U n=0 

00 



W-\e) = J2(-c) m X m (6.14) 



m=0 



the general term of the series is of the form (— e) n+m R n {a s {t) / u) and is reduced to a power 
series in a s (t) by repeated application of the identity 



-^ t )p(^) 9 « (i) = Q ,( iho (^)^) ff « (i ) _ A ^ {t)p ^ g m m 

(6". 15) 

valid for any function p, where 

1 + + (616) 

Eq. ()6.15|) means that multiplication by — e corresponds to the operator 7o(a s (£)/u;) 
after integration by parts. It follows that the quark anomalous dimension is given by 



n=0 



(6.17) 

which is the resummed expression we were looking for. 

The subtlety of the result (|6.17|) is that, in order to get a resummed formula in a s (t)/u) of 
the e-independent 7 qg , we need an all order e-expansion of the fc-factorisation formula (j6.12j) . 
The fact that terms of order e n generate finite contributions to both coefficient and anomalous 
dimension is somewhat similar to what already noticed in the gluon channel, and is typical 
of the minimal subtraction recipe. However, while higher orders in e correspond to higher 
subleading logl/x resummation levels in the gluon case, here they just correspond to higher 
orders in the a s (t) /^-expansion of 7qg . 

In order to use Eq. ()6.17|) we need to understand the action of the 7o (l + -D) -1 operator. A 
simple example is provided by setting R n = 5 n0 and 7o = a B (t)/u. By noting that D[a s (t)] n = 
n[a s (t)] n , we obtain 



(MS) «b(*) TT^Mn f a B (t)\ f A «rn*( a ^( t )\ a s(t)T R ( ,M1 1|M 



n=Q 



(6.18) 



where Kb denotes the Borel transform of !K rat , which is simply obtained from Eq. (J6.11|) . As 
noticed in [6], Eq. (j6.18|) provides all rational coefficients occurring in the resummed formula 

r (MS) 

for7 qg ; . 

It is also possible to provide an expression for (jQ7j) involving only quadratures, provided 
the functions R n are given to all orders. In App. [D] the result is derived in terms of the 
intermediate function 
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where Rb is the Borel transform of K tva,n (e)R(a s /uj } e) in the — e variable. The final result is 
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(6.20) 



where we notice the exponentials of 1)6.18)1 occurring in a more general framework. 



7 Discussion 

The main results of this paper are the renormalisation group (RG) factorisation of the BFKL 
equation in 4 + 2e dimensions at NLx level, and the relation of the Qo-scheme to the MS- 
scheme at the same level of accuracy. The collinear factorisation has been proved by solving 
the BFKL equation with a cutoff Qq = /x 2 e _T , and by showing that — in the on-shell limit 
for the initial gluon (Qq = 0) — the solution admits a RG representation with a well-defined 
anomalous dimension %(a B ;uj), as given by Eqs. (|3.10jh ()4.5|) and ()5.10|) at various levels of 
accuracy. This representation shows exponentiated IR poles ~ ^ (^r) > which are factorised 
in the IR-free regime with a s (t) < e/b and a s (— 00) = 0. Then, by the e-expansion of 7 e (o; s ; u) 
in Eqs. ()4.3|) and f|5.9|) . we are able to define the minimal subtraction scheme, to switch to the 
UV-free regime, and to find the transformation factor R(a s ; u) at NLx level. 

We discover in this way that the coefficient R is due to the product of a fluctuation factor 
N, which can be calculated at e = 0, and of a dynamical factor 9^ = RfN which reflects 
the ^-dependence of 7 e (o: s ;u;), expanded around e = ba s . Furthermore, the e-dependent 7 
directly provides subleading contributions to the MS-scheme anomalous dimension, encoded in 
7 £ =te s (a s ; w) (Eqs. dHHD, (IHTTf and (fHH2D )- In this way, the difference 7 (M5 > - 7W ) between 
the two schemes is here calculated up to NNLx level. Therefore, the e-dependence of the kernel 
is transmuted into a subleading a s -dependence for the MS anomalous dimension. 

A similar transmutation phenomenon occurs in the case of quark-gluon mixing. In addition, 
due to the different form of the scheme- changing transformation ()6.2|) . the e-dependence induces 

leading a s /d;-dependence into 7qg^ — which has to be disentangled from the pro cess- dependent 
contributions to Cq^ . Therefore, the knowledge of the e-dependence becomes increasingly 
important for the full singlet evolution. 

The above results are not directly applicable to the doubly resummed approach [14] nor, as 
far as we understand, to that of [15]. However, they provide some hints towards an improved 
scheme-changing transformation. First of all, the analysis of Sec. E] and App. [B] shows that 
there is a universal part in the e-dependence of K^ (7) (the p-scheme defining kernel) which is 
encoded in the collinear pole at 7 = — e, which in turn comes from an off-shell generalisation 
of the well-known P qg splitting function. It is then conceivable that a similar analysis can be 



performed for resummed models as well, in order to control the leading part of the scheme- 
changing transformation and its mixing. 

Furthermore, since the higher order e-dependence affects R at subleading level, it is con- 
ceivable that most of the normalisation change comes from the already known (or generalised) 
leading part. A preliminary analysis in this direction is under way [20]. Hopefully, this will 
soon lead to a realistic comparison of the resummed approach to experimental data. 
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Appendices 

A Solutions of the difference equations in 7-space 
A.l Asymptotic solution 

In Sec. |2] we showed that the solution of both the homogeneous and inhomogeneous equa- 
tion ()2.8|) in 7-space are determined by difference equations whose structure is, in both cases, 
given by (cf. Eqs. (I2~T5J) and (jQZj) ) 

S £ ( 1 + e)-S £ ( 1 )=L £ ( 1 ). (A.l) 

Evidently, Eq. (jA.l|) determines S £ {^f) up to an additive periodic function of 7 of period e. Here 
we want to provide a solution for S e in terms of the known function L £ . Due to the non-local 
form in 7 of Eq. (jA.lJ) . we expect S e to be determined by an infinite number of 7-derivatives 
of L e . In fact, <9 7 is the generator of 7-translations, and exp(e<9 7 ) the translation operator: 

[exp(£<3 7 )S e ](7) = S £ (7 + s), (A.2) 
so that, by substituting Eq. (jA.2|) into Eq. (jA.l|) we obtain the formal solution 

S £ = [exp(e9 7 )-l] _1 L £ . (A.3) 



Since we are not able to provide the resolvent in Eq. (jA.3|) in closed form, we look for a 
perturbative solution in e of S £ . Given that exp(e<9 7 ) — 1 ~ e<9 7 for e — > 0, we can rewrite 
the r.h.s. of Eq. (jA.3|) in terms of an operator close to the identity by applying ed-y to the left, 
obtaining 

C £ ee ed,S £ = £0 \ L E = £ ^ n L[ n) (7) , (A.4) 

v " n=0 

where we have adopted the definition (j2.12j) and used the generating function (|2.15jl of the 
Bernoulli numbers. This proves Eq. (j2.14j) . 

A remark concerning the character of the Series in Eq. (|A.4|) : since the asymptotic behaviour 
of the Bernoulli numbers is B n ~ n\, it turns out that the series has a non- vanishing radius 
of convergence if and only if the series Yl n L £ {^)x n /n\ has infinite radius of convergence, i.e., 
L £ ( , y) is holomorphic in the whole 7-plane. If L e {^y) has at least one singularity at finite values 
of 7, then the series in Eq. f|A.4|) is an asymptotic series of C e (j) for e — > 0. 



A straightforward example is provided by the Euler ip function which satisfies ip(z + 1) — 
ij}(z) = 1/z. By rescaling z = 7/e and = S r e (7), we reproduce Eq. (|A.1|) with L e (7) = 
e/7, hence Lg (7) = e(— l)"n!/7 n=1 which, substituted in Eq. ()A.4|) yields = £e(t) = 

J2n=o B n (-l) n (e/ 1 ) n+1 and finally ij)(z) = \ogz + J2n=i B n (-l) n+1 /(nz n ) + c, which is exactly 
the asymptotic expansion of the if) function for large \z\, provided we set the constant of 
integration c = 0. Note that also the function z i— > — z) satisfies the same difference 
equation as i[>(z), an d i n f ac t they differ by the periodic function 7rcotg(7rz). The asymptotic 
expansion of ip{l — z) differs from that of if)(z) by the constant — i7rsign($5(,2)) and is therefore 
still consistent with our series representation. In practice, this kind of ambiguities requires 
the choice of a Riemann sheet for the asymptotic solution S^j), which has one or several 
branch-points. This is done on the basis of the regularity requirements for the BFKL solution 
in the UV or IR region. For instance, the saddle point (|2.35J1 is dominant in the Riemann sheet 
relevant for the UV regular solution of type (|2.34j) . 



A. 2 Iterative solution 

The BFKL equation can be solved through the iterative method, both in t- and 7-space. In 
this section we want to prove the equivalence of these two procedures. 

Let us first review the iterative method in t-space. From the integral equation (J2.2j) for 
J-'eik) we derive the corresponding equation for -F £ (t) = .F e (fe): 



J dt'K e (t- 0^(0 



(A.5) 



u 

where the rescaled and azimuthally averaged kernel 

KJt - f) = -A^- (XJk, k')) = f- / d 1+2e k X e (k, k') (A.6) 

is a dimensionless scale- invariant kernel whose eigenf unctions are exponentials e 7 * : < 3?7 < 1, 
and the corresponding eigenvalue function is 

Xe{l) = e~ 7 * J dt' K £ {t - t)&* , (0 < 3fry < 1) . (A.7) 

The iterative solution of Eq. (jA.5|) is then given by 

00 

^) = E^ 1] (*), (A.8) 

n=0 

where 

FM(t) = ^e £t , (A.9a) 

^N( t ) = fV J dt , Ke(t _ t ^[n-i] (t ') = (% £ *) n+1 f[ xe(ke) (A.9b) 



k=l 



so as to yield Eq. (J2.4j) . Note however that the action of the kernel K s on Tf" ^ ~ e n£t in 
Eq. (JA.9bj) is defined only if ne < 1. For n > 1/e the functions TT ^ do not belong anymore 
to the domain of K £ and the iterative procedure breaks up. Nevertheless, from a perturbative 
point of view, Eq. (|2.4jl is meaningful, because for any given order n the coefficient of a s (t) n+1 
is well-defined and analytic in a non-empty open interval ]0, l/n[ of e- values. 



In 7-space, on the other hand, the iterative solution (|2.23|) has been obtained by using an 
initial condition with an infrared cutoff T — see Eqs. (I2.9I2.10JI in comparison with Eq. ()A.9aJ) 
- so as to deal with analytic (non symbolic) functions. 6 In order to prove the equivalence of 
Eq. (J2.4J) and Eq. (J2.23J) . we rewrite the latter in the form 

00 

/.ft) =£/.%), (A. 10) 

n=0 

/a. \ n+1 e [7 - (n+1)£]T A 

^w=(2) ^-rcn^-fe), (A.11) 



k=l 



where we have used the explicit form (J2.9J) of In performing the inverse Fourier trans- 
form (J2.7J) . the real part K7 = c of the integration path must lie to the right of the UV 
singularities of /J™' (7) at 7 = me : m = 1, ■■■ ,n + 1 and to the left of the IR ones at 
7 = 1 + me : m = 1, • • • ,n + 1. Therefore, also in 7-space the iterative solution appears 
to be defined only for n < 1/e. Here, however, we can deform the integration path and extend 
the validity of the 7-representation to higher values of n, though this requires crossing the real 
axis at larger and larger values of 7 > ne and to properly exclude the IR poles of the eigenvalue 
functions Xe{l ~ for 7 — m £ > 1- 

We thus calculate, for n < 1/e and (1 + n)e < c < 1 + e 

m*) = / ^ = ^ +i, - T I ^ — -7—7-pr" n x.(7 - *«) . 

J^=c 2vri V u J J^ 1=c 2m 7 - (n + l)e ^ 

(A.12) 

For t > —T we can close the contour path to the left, and pick up the rightmost pole of at 
7 = (n+l)e and other poles of the Xe's at 7 = me : m < n + 1. Each pole contributes with a t— 
and T-dependent exponential of the form exp[— (n + 1 — m)eT + met] : m < n + 1. Therefore, 
the residue of the rightmost pole (m = n + 1) is T-independent, while the other residues are 
suppressed as integer powers of exp(— eT) and vanish in the T — > +00 limit. In conclusion 

(- /a \ n+1 n 
— ) n^((^ + l-^) + ^(e- £T ) (A.13) 
w ^ fc=i 

and, in the T — ► +00 limit we reproduce the t-space iteration (|A.9b|l . 



A. 3 Calculation of the p-function 

The calculation of the function p e (j) obeying Eq. (J2.27J) requires a special care with respect 
to the calculation of ^(7) in Eq. (J2.13J) . even at leading e accuracy, because the r.h.s. of 
Eq. (J2.13J) is finite for e — >• 0, while that of Eq. (J2.27J) has an essential singularity. By using the 
expression (J2.12J) for h £ {j) into the r.h.s. of Eq. (J2.27JI . according to Eq. (|2.14J) we expect 



B 



^(7) = Eif^r )(7)e " L£(7)e 



-^(7) p -U 7 ^(V)dV' 



(A.14) 



n=0 



In computing the logarithmic 7-derivative of the product in square brackets, the large contri- 
butions are represented by T (from the exponent of f^°\ cf. Eq. (J2.9JI ). and by the derivative 

6 The Fourier transform of an exponential is a S function. 



of the "action", which is of O (1/e). The other contributions are T- independent and regular in 
e, so that for higher derivatives one gets 



d" 



/ (o) (7)e -^(7) e -ir£ £ (7Wj =/ (o) (7)e -L £ ( 7)e -ir£ £( y)dy (V_I £e(T) y x [i + o(e)] . 

(A.15) 



e 



Eq. (JA.14j) becomes then 



oo 



^ 7 p £ ( 7 ) = /W( 7 )e-^We-*W £ ^ [eT - L ( 7 )] n x [l + (e) ] 



n=0 



= /(»)( 7 )e-M7)e-4(7)^_^ x [l + (e) ] , (A.16) 

where in the last equality use has been made of Eq. (J2.15|) . Finally, after integration in 7, one 
ends up with Eq. 

B Characteristic functions of the quark kernels 

In this section we want to compute the characteristic functions of the g — > q kernels used 
to define the MS and Qo quark densities in the high-energy regime according to Eq. (|6.1|) . 
These kernels are naturally provided in (x, fe)-space, and are related to their cu-moments in 
Eqs. (|6.U6.3j) through a standard Mellin transform. In the collinear limit, they turn out to be 
all determined by the same off-shell generalisation of the P qg splitting function. 

B.l The Catani-Hautmann kernel ij( CH ) 

The g — > q kernel used to define the quark density in the high-energy regime is given in 
Eq. (4.8,4.9,0.6) of Ref. [6] (we call f their variable z): 



2, 1 , 1,2 s 

(B.l) 

Ref. [6] 



^ CH, (4KM4 



where t = logQ 2 //i 2 along this section. We want to compute the double Mellin transform (cf. 
Eq. (Q) 

^^* tf )f d «-/f (S) 7fffH, (4) ■ (a2) 

The result at u = can be read from Eq. (C.9) of Ref. [6] and in our notations is 

«f (7) - *?»><M> - + - 3, + ^ T + :%:^.7> ■ (B.3) 
The coefficient of the collinear pole at 7 = — e is 

v; e k J 2tt 3 (l + 2e)(l + |£) r(l + 2e) v ; 



(T) 

Figure 2: The process 7*g — > qqg used for computing "K\ . 
B.2 The transverse kernel H( T ' 

We define the Qq quark density as the high-energy limit of the quark density in the DIS- 
scheme, i.e., q( Qo ^ = F 2 \ high energy, the corresponding process being depicted in Fig. |21 We 
choose a reference frame in which the incoming momenta q and p define the longitudinal plane. 
Then, in the high energy limit, the momenta entering the blob can be decomposed as follows: 



q = n — xp ; k = zp + k . 



(B.5) 



where n and p are light-like momenta spanning the longitudinal plane orthogonal to the trans- 
verse plane. In the latter we use the euclidean representation fc^fe^ = — k 2 < 0. The remaining 
momenta have the following Sudakov decomposition: 



I = p q — q = k — Pq = —f3n + ap + I 
Pq — q + I — (1 — (3)n + (a — x)p + I , 
Pq = k — I = (3n + (z — ot)p + k — I . 



(B.6a) 
(B.6b) 
(B.6c) 



We then decompose the structure function F 2 in transverse and longitudinal part, each of which 
is given the factorisation formula 1)6.1)1 with an appropriate kernel 



a. 



9 2 T R Q 2 



(2n] 
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3+2e 



£k jo 



dp 



l2+2s 
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I 2 



(k-l) 

[i-p)q 2 m 2 



Ci{/3,Q,l,k) 



where ^ = x/z and 



C L = 8P 2 (1-P) 2 Q' 



Ct = Pq,(P) 



I 



[D(l) 
I - k 



2 



D(l) D(l-k)_ 



D(l) = l 2 + P(l-P)Q 2 ; PqJ(3) 



I) 



(3 2 + (l-(3) 2 + e 



1 + e 

The Mellin transform are again computed according to the definition ()B.2 
integration in £ is easily done by exploiting the 5-function: 

7-1 



« s (t)M«( W , 7 ) 



g 2 T Rl (j + e) 



(27T) 



3+2e 



d/5d 



2+2e 



I 



dk' 



Q 2 



i 



i z 

(i-P)Q 2 



k-l) 



(B.7) 

= L,T), 

(B.8) 
(B.9) 

The first 



a. 

(B.10) 



The case e = has been computed exactly in Ref. [19]: 

s4V 7) = 9 J* 7 ' S ^ N = w ' 7)l ^- m • (R11) 

Z7T u; em tl! s 

At £ 7^ we are mainly interested in the u = expression of the transverse kernel, which 
provides the singular anomalous dimension behaviour 



a, 



.mPi^^^mPion) (B.i2) 



9 2 Tnl{l + e) f 1 AQ D _ /" dfc 2 /fc 2 ^" 1 



dj9 P q7 (/?) / ^ ( ^ ) / d 2+2£ Z 



I l-k 



1 2 



D(l) D{l-k) 



The Z integral can be performed by expanding the square, shifting integration variable l — k*->l 
in the square of the second term inside the brackets, collecting all terms into a single fraction 
and using Feynman's parametrisation of the common denominator 

i _ yi 2(1 -A) 



D 2 (l)D(l - fc) 



/ dA < 2, M2Y3 > P = l~^k, M 2 = (3(1 - (3)Q 2 + \(1 - \)k 2 . 
Jo {P + M ) 

(B.13) 



The corresponding numerator expressed in terms of the new variables reads 

(1 - X)k 2 p 2 - 2A(fc • p) 2 + Afc 2 M 2 + odd terms in p . (B.14) 

Then, by means of standard integrals in 2 + 2e dimensions, one can perform the p integrals 
obtaining additional (3- A- and fc 2 -dependent factors. One can now perform the k 2 integration 

d*V*"Y- e - < y B(7 , 1 - T - e) W-^)ir 1 , (B.15) 



k 2 \Q 2 J M 2 (!- £ ) ^ v " ' ' [A(1-A)]t 

which provides, among other things, the e-running of the coupling. Then we evaluate the A 
integral 

1 dA (1 - A)[(l + e)(l - A) - £ A][A(1 - A)] -7 = f±lZ2 (B.16) 

and the (3 integral 

q7VPnPV r(2 + 2 7 + 2e) (i +e )( 7 + e ) v 7 

obtaining the final result 

r K{ T )( \ = Tr e^)r(i + 7 )r 2 (i - 7 )r 2 (i + 7 + e)r(i - 7 - g) 

£ l7j Tr r(2 - 2 7 )r(2 + 2 7 + 2e) 

x (2-7 + ^(1 + ^) + 1 + ^ + ^1 (R18) 
(l+e)(3-2 7 ) ' 

The coefficient of the collinear pole of %£ ( 7 ) at 7 = — e equals twice the analogous quantity 
%(e) of Eq. (jB.4|) obtained from the Catani-Hautmann kernel, because here we consider both 
the quark and antiquark contributions in the fermion loop, while in the kernel K qg given in (jB.l|) 
only the quark (or antiquark) contributes. 



B.3 Universal collinear behaviour and off-shell splitting function 

The fact that the residue at 7 = — e of !Ke ( r y)/ r y('y + s) is independent of p is a simple 
consequence of the universal collinear behaviour of in the region I 2 ~ k 2 <C Q 2 , when 
either I or k — I are collinear with k. In this region we can solve the 5 function in Eq. (jB.7|) for 
small values of (3, namely 



«1, 



(B.19) 



i-e <2 2 

or instead in the symmetrical region of 1 — (3 small, with I replaced by k — l. Taking the region 



f|B.19j) . which corresponds to a collinear quark, we see that the longitudinal part is power 
suppressed, so that it plays no role in the following, and does not contribute to the 7 = — e 
pole. Furthermore, we can set P q7 (/3) = 1 in the transverse part, with the denominators 



(1 - i)D{l) ~ (1 - Ol 2 + £(* -l) 2 = f + ^- Ok 2 = D 
(l-£)£>(Z-fc)~(fc-Z) 2 , (l = l-ik). 

By replacing (lB~T9ll and (lB~20l) in (IB~7|) we obtain 



9 2 T, 



(27r) 3 + 2e fc 2 



d 2+2£ l 



~ l 2 (k — I) 2 + 21 ■ (k — l)D + D 2 



D 2 



(B.20a) 
(B.20b) 



(B.21) 



Introducing the boost invariant transverse momentum I = I — £k, the numerator in Eq. ()B.21|) 
can be recast in the form 



l2 



Vk 2 - 4^(1 - fl(fc • l ? + 4 ^(1 - 2 (fc 2 ) 2 + 4^(1 - 0(1 - 2i)k 2 {k ■ I) 



(B.22) 



~2 



which, after azimuthal averaging in I at fixed I and k , yields (7 

„ 2 „Uf\r>2 ,?2 / ~2 \ e 



2 

max 



/ (1 2 \ e £ ^(i) rH£)Q 2 A] 2 /f\ , r.2 

^k 2 J-2nT(l + e)J f \&) r <*\?>?> 



= A(0Q 2 ) 

e), (Z 2 ,(fc-Z) 2 «Q 2 ) 

r(CH) 



(B.23) 



This behaviour is identical to that of the Catani-Hautmann kernel H £ , with the transverse 
momentum dependent splitting function 



p<f>)(£ 
1 qg v^' 



k 2 ,e) 



1 



e + (l-Q 2 + e 
l + e 



+ 4f(l-0V 



K 



(B.24) 

which thereby acquires a universal meaning. The corresponding contribution to < K ( f\uj^) 
becomes 



7(7 + e) 



coll 



d^ e 

e eV(l) 



dfc 2 /fe 



„2\ 7 



27rr(l + £) 



de e 



^ 2 df // 2x7+£ 



(B.25) 
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7~-e 
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3 e^(i) 



1 + e 27rr(l + e) J 



de c 
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"dK 2 



d K 2 
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(B.26) 
(B.27) 



7=-e 



We see that the residue at the 7 = — e pole is independent of the details of the I phase space 
I < \{i)Q 2 and is given by a simple (analytically continued) moment of the P q g splitting 
function. We obtain, for u = 0, the explicit expression 



c 



JC(e) = — — / d£ / dn 2 (« 



2ttT(1 + s 



jo 



-2 t..2\-e 



2 m^*'> K ' e ^ 



8k 



(B.28) 



t r 2 i_ + e e^)r 2 (i + g )r(i- £ ) 

2tt 3 (l + 2e)(l + |e) r(l + 2e) 1 ' J 

The above analysis shows that, by combining fc-factorisation with the collinear behaviour it is 
possible to define an off-shell splitting function which is process-independent. The residue at 
the pole 7 = — e is a particular moment of such splitting function. 

C O (^-corrections to the gluon density 

In order to compute the corrections of relative order ba s to 9t and to 

7 (MS) we 

have to calculate 

the corrections of order O (e) — i.e., of relative order O (e 2 ) — to the exponent of the solution 
(PEE} of the (unintegrated) gluon density T £ . 

C.l Calculation of the saddle point fluctuations 

These corrections can be computed by means of a saddle point expansion of the 7-representation 

F £ (t) = M £ J d 7 e^ +& W. (C.l) 
Here the "action" S e has to be expanded to O (e) as done in Eq. (j2.17|) . and the normalisation 



can be inferred from Eq. (|2.34|) . which also sets to the lower bound of the integral in the 
action. We split the exponent of Eq. (|C.1|) into an e-regular part 



plus a leading-e: phase 



EM = ~LM + ±L'M + O (e 2 ) (C.3) 

E £ ( 1 )^ 1 t + - f d 7 'M7'), (C4) 
£ Jo 

whose stationarity condition determines the saddle point 7 = % (cf. Eq. ()4.2j) ). and the size of 
the fluctuations 

a ' " = ■ (c ' 5) 

By expanding around the saddle point and factoring out of the integral the value of the integrand 
at the saddle point, we can write 



T £ (t) = M e e^ &(%) V27r<7 7 ( exp 



,k=3 k=l 



, (C6) 



where the "average" (• • •) is defined by 



</(7)> 



2na^ 



exp 



(7 



2<r* 



(C.7) 



In Eq. (jC6j) the sum of the E^ starts from third order because the zeroth order has been 
taken out of the integral, the first order is zero by definition of the saddle point condition, and 
the second order provides the gaussian measure in the "average integral" (|C.7[) . 

The product out of the "average" in Eq. (|C.6|) — by using the expressions APSED for M £ 
and a 7 and Eqs. (j2.17p3.7|4.6j) for the exponential — yields 



2lXO^ 



1 



exp 



\/-Xe(le) 

i.e., the gluon density of Eq. (j3.9|) times a correction factor 



dr 7 e (r) 



(C.8) 



ensile) 



l+e-L' £ (% 



O (e 2 ) . 



(C.9) 



The remaining corrections are found from the "average integral" . Note first that the "average" 
of powers of 7 — 7 e vanishes for odd powers, and is given by 



<(7-7 £ ) 2n ) = (2n-l)!! 



a 



2n 



{2n- 1)!! 

[-L'MT 



(CIO) 



for even powers, thus providing higher and higher orders in e as n increases. However, the coeffi- 
cients of the saddle point expansion present ^-singularities, due to the leading phase derivatives 
E^ ~ 1/e. When one expands the exponential inside the "average", the m-th order (which 
is of 0(l/e m )) contains powers (7 — %Y only for j > 3m. Therefore the power in e of the 
fluctuations (jC.lOj) increases faster than the inverse power in e of the order of expansion of the 
exponential. In practice, if one is interested in computing the "average" to, say, l-th order in 
e, it suffices to expand the exponential up to order m < 21 and to take terms in the sum up to 
k < 21 + 3 — m. In the present case, we are interested to the first e-corrections, namely / = 1, 
and we need therefore the first order expansion of the exponential up to k = 4 and the second 
order up to k = 3. By denoting 5 = (7 — %) we have 

/ R(2) _ R(4) . 



1 + ^r 6 + ^r 6 + 2 



= l + e 
where we have used 
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(C.12) 



and the dependence on e and 7, has been omitted. Finally, the full O (e) correction to the 
action ()4.8|) is obtained by multiplying the saddle point fluctuations of Eq. (jC.llJ) by the first 
order correction to the action at the saddle point in Eq. ()C.9j) . 



C.2 Cancellation of the singularities 

The terms of the O (e)-correction to the gluon density Si present various singularities deriving 
from those of the function L' £b {j) written in Eq. (j4.11|) : (1) poles at % = 0, i.e., at a s (t) = 0; 



(2) poles at Xe(le) = bu/e, i.e., at a a (t) = e/b; (3) poles at e = because of the 6a s /e-expansion 



in the denominator of L 
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In order to better understand the role of these singularities, we rewrite Si in terms of the 



inverse of L' 
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and of the eigenvalue functions Xe and its derivatives. By computing the derivatives of A 

XL"? 



(C.13) 



L" 

W) 



X'- 

we deduce 

(L'T 
{L'f 

so as to obtain 



1 - —A , A" 



L" 



{L'f (L>y 



Xe i 
'X'e 



Xe 
X'e 



A (C.14) 



X 



- -2^f + A ^ 



Xe 



2 y" 
T-3^ + A 



//\ 2 



5l 



24 X' 



i ^xjA' 

X' 



12 Vx £ 



A 



(C.15) 



(C.16) 



which agrees with Eq. ()4.12|) after replacing the expression (|C.13|) of A and evaluating at 7 = j £ . 



D Expression of 7 qg by quadratures 

In this appendix we want to show how it is possible to give an expression for the NLx resum- 
mation of 7qj^ IS ' ) in terms of the universal functions R(a s /u,e) and "H(s) up to quadratures. 

We start from Eq. (j6.17|) and notice that, besides the overall factor a s (t), all the remaining 
dependence on a 8 (t) occurs through the ratio a = a(t) = a B (t)/uj. In fact, also the operator D 
defined in Eq. f!6.16|) can be written as D = ad a . Let us then introduce the operator 

6 = 7o [l + D]- 1 = l0 (a)[d a o a]' 1 (D.l) 

and, by using the expansions (jbM3lbM4|) . rewrite Eq. (|(j.!7j) in the form 
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The first step is to derive an expression for arbitrary powers of O. We solve this problem 
by considering the resolvent of O 

00 

Z (5 (A)EE[l-AO]- 1 = ^A fc O fc , (D.3) 

fc=0 

where A is a complex parameter. The resolvent is the generating function of the powers of O, 
in the sense that 

1 d n 

n\ dA n ov ' 

An expression for Zq(X) is obtained by solving the equation 

y(a) = [Z d (X)x](a) <^ x = [1 - XO}y = y - A 7 q[1 + D^'y (D.5) 



(D.4) 

A=0 



where y(a) is the unknown function to be determined in terms of x(a). By introducing the 
function w(a) such that y = [1 + D]w = d a (aw), we obtain a first order differential equation 
for the function aw (a) which can be easily solved. The general solution of Eq. (|D.5|) is 



[Zq{\)x\{q) = d a da' exp 
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The lower bound of the a'-integral is 0. It can be determined by comparing the expression of 
O we obtain from Eqs. (|D.4|) and (|D.6|) 
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with the integral representation of the operator O that can be obtained by the series expansion 
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The second step is to perform the sum in n in Eq. (jD.2|) . By using again Eqs. (jD.4|) and 
ID. 61) we obtain 
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where we have noticed that the sum inside the integral is just the Borel transform of the function 
!K tran (e)i?(a, e) in the (-^-variable. 

The third step is to perform the sum in m in Eq. (jD.2|) . Here we exploit the particular 
structure of the coefficients IK m of being a linear combination of powers, since CK rat (e) (cf. 
Eq. (JbMl)0 is a linear combination of simple fractions: 
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so that, by using the second equality of Eq. (jD.3|) . we obtain just a linear combination of 
resolvents with parameters A = 2 and A = 2/3: 
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Finally, by inserting Eqs. (lP~9|l and (IP.llj) in Eq. (lP~2|) yields 
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